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We characterize the conditions under which a translationaUy invariant matrix product state (MPS) 
is invariant under local transformations. This allows us to relate the symmetry group of a given 
state to the symmetry group of a simple tensor. We exploit this result in order to prove and extend 
a version of the Lieb-Schultz-Mattis theorem, one of the basic results in many-body physics, in the 
context of MPS. We illustrate the results with an exhaustive search of 5[/(2)-invariant two-body 
Hamiltonians which have such MPS as exact ground states or excitations. 



I. INTRODUCTION 

Matrix Product States (MPS) [D I2 encapsulate many 
of the physical properties of quantum spin chains. Of par- 
ticular interest in various physical contexts is the subset 
of translationaUy invariant (TI) MPS, originally intro- 
duced as finitely correlated states [Tl. Their importance 
stems from the fact that with a simple tensor. A, one can 
fully describe relevant states of N spins, which, at least 
in principle, should require to deal with an exponential 
number of parameters when written in a basis in the cor- 
responding Hilbert space TC^^ . Thus, all the physical 
properties of such states are contained in A. It is there- 
fore important to obtain methods to extract the physical 
properties directly from such a tensor, without having to 
resort to H®'^ . 

An important physical property of a TI state, vp, is the 
symmetry group under which it is invariant. That is, the 
group G such that 
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where g (z G and Ug is a unitary representation on Ti.. In 
a recent paper we showed that for certain kind of MPS 
(those fulfilling the so-called injectivity condition [TJ[2]), 
this symmetry group is uniquely determined by the sym- 
metry group of A (with a tensor product representation) . 
Roughly speaking this means that by studying the sym- 
metries of A we can obtain those for the whole state 
This result allows us, for example, to shed a new per- 
spective into string order a key concept in strongly 
correlated states in many-body quantum systems. 

Another relevant property of MPS is that they are all 
exact ground states of short-range interacting (frustra- 
tion free) Hamiltonians [11 [2]. In particular, for every 
TIMPS we can always build a (so-called) 'parent' Hamil- 
tonian for which it is the ground state. Of particular in- 
terests are TIMPS with two-body parent Hamiltonians; 
that is, whose parent Hamiltonian consist of two-body 
interactions only. And among those, the ones which 
have a large symmetry group, like SU{2). The reason 
is that those are the ones that naturally appear in con- 
densed matter problems. Two prominent examples are 
the AKLT |4] and the Majumdar-Gosh [5] states, who 
have two-body parent Hamiltonians with SU(2) symme- 



try. They have served as toy models to understand cer- 
tain physical behavior in real physical systems, like the 
existence of a Haldane gap [6] in spin chains with integer 
spin, or the phenomenon of dimerization [51, respectively. 
Despite their key role in the understanding of spin chains, 
there are very few other examples known of TIMPS with 
SU(2) symmetry and with a two-body parent Hamilto- 
nian [niHjin]. 

In this work we first generalize the results of Ref. [3] to 
arbitrary TIMPS. This enables us to derive some generic 
properties about those states, as well as to obtain a sim- 
ple proof for a version of the Lieb-Schultz-Mattis theorem 
[TU]. This celebrated theorem states that all Hamiltoni- 
ans with SU{2) symmetry are gapless for semi-integer 
spin (dim(7i) = n -I- 1/2, n = 0, 1, . . .). In our case, we 
can prove that all TIMPS corresponding to systems with 
semi-integer spins cannot be the unique ground state of a 
local frustration-free Hamiltonian. Furthermore, we can 
extend the proof to other groups, like U{1) for spin 1/2 
systems, and find counterexamples for this last case when 
the spin is 5/2 or larger. 

In the second part of our work we concentrate on MPS 
that are eigenstates (not necessarily grounds states) of a 
(so-called 'parent') Hamiltonian which has SU{2) sym- 
metry and contains two-body interactions only. We find 
other families of Hamiltonians beyond the well-known 
AKLT and Majumdar-Gosh with those features. Fur- 
thermore, we find the first examples of MPS that cor- 
respond to excited states of SU (2)-invariant Hamiltoni- 
ans. There is a new example of state with spin 1, which 
is never the ground state of any frustration free SU (2)- 
invariant two-body hamiltonian. In order to make a sys- 
tematic search of all those MPS we develop a simple tech- 
nique that allows for a numerical systematic search. 

This paper is organized as follows. In Section [H] we re- 
view some of the basic properties of TIMPS and establish 
the notation that will be needed in the following. In Sec- 
tion HI we establish the relation between the symmetry 



group of a TIMPS and that of the tensor A defining the 
MPS. For continuous symmetries, such as SU{2), we will 
see that the set of symmetric TIMPS is intimately related 
to the set of Clebsch-Gordan coefficients. Section lTvl then 
provides an MPS version of the Lieb-Schultz-Mattis the- 
orem and in Section |V] we give a detailed investigation of 
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SU{2) symmetric TIMPS which are eigenstates of two- 
body Hamiltonians. 

II. MATRIX PRODUCT STATES 

Let us consider a system with periodic boundary con- 
ditions of N (large but finite) sites, each of them with 
an associated d-dimensional Hilbert space. A transla- 
tionally invariant MPS on this system can be defined 
with a valence bond construction in the following way: 
Let us consider another couple of D dimensional ancil- 
lary/virtual Hilbert spaces associated to each site and 
connected to the real/physical d dimensional space by a 
map A — J2i a ^i,a/3 N) {o^Pl ■ Then, by introducing max- 
imally entangled states connecting every pair of neigh- 
boring virtual Hilbert spaces (usually called entangled 
bonds), it is not difficult to prove that the state can be 
written as 

1$)== tT[A,,---A,^]\ii---iN) 

ii,...,ijv 

where we call the matrices K, — {Ai E Aijj , i ^ 1, . . . ,d} 
Kraus operators. A way to work simultaneously with all 
of them is to define the map 

F^^A,®|z). (2) 

i 




FIG. 1: This figure represents the MPS construction. A pair of 
virtual spins which are connected to their neighbors via a maxi- 
mally entangled state |f! >= \oi ot > are mapped into the 
physical spins (below). All properties of the state originate from 
the mapping between physical and virtual system. 

For each MPS there exists a canonical form [2J The- 
orem in. 7, Lemma IV.4] which assures that one may 
choose all matrices Ai with a block diagonal structure 
|22j . in such a way that after gathering enough spins to- 
gether, the Kraus operators fulfil: 



Property 1 {Span property). The set of products V = 
{Ai-^ ■ ■ ■ Ai^}, with n the collected spins, spans the vec- 
tor space of all matrices with the same block diagonal 
structure. 

It is an open conjecture stated in and verified in 
many particular cases, that an upper bound for the num- 
ber of sites which have to be gathered to achieve property 
[T] depends only on the dimension D of the Kraus opera- 
tors. When there is only one block in the above canonical 
decomposition the MPS is usually called injective, since 
the linear operator mapping boundary conditions to the 
resulting states is indeed injective [TJ|5] when taking suf- 
ficiently many particles. The definition reads: 

Property 2 (Injectivity). There exists n such that the 
mapr„(X) = Eii,...,j„ ^r{XAi^ ■ ■■A,J\ii ■■■in) is injec- 
tive. 

For each MPS {ip) one can construct a Hamiltonian, 
called parent Hamiltonian, for which ji/i) is an eigenstate 
with eigenvalue . 

Definition 3 {Parent Hamiltonian). Let p^''^ be the re- 
duced density matrix of |?/;) for k particles (k will be 
called the interaction length of the parent Hamiltonian). 
Let us suppose that {|fi)}r=i, with r > 1, is an or- 
thonormal basis for ker [p'*^'] . Taking any linear com- 
bination of projectors h(d) = X)I=i define 
H = J2i Ti{h)®lrcst, where t; is the translation operator. 

If ai > 0, then the Hamiltonian is positive semidefi- 
nite and \ip) is indeed a ground state. Moreover H is 
frustration free, since \tp) minimizes the energy locally. 
Injectivity has now a deep physical significance. If it is 
reached for n particles and every Oi > 0, it ensures that 
the MPS is the only ground state of its (n+l)-local parent 
Hamiltonian, that it is an exponentially clustering state 
and that there is a gap above the ground state energy 

WE]- . 

In this work we will focus on symmetries of states in- 
stead of Hamiltonians. There is however a close connec- 
tion between the two approaches. On the one hand, it is 
clear that the unique ground state of a symmetric Hamil- 
tonian has to keep the symmetry. On the other hand, we 
have the following 

Proposition 4. If an MPS \ip) is invariant under a rep- 
resentation of a group, one can choose its parent Hamil- 
tonian H invariant under the same representation. 

To see that it is enough to notice that the symmetry in 
the state ([T| implies the invariance of ker [p^*^-*] under the 
same symmetry. Symmetrizing ker [p^*^-*] (i.e., averaging 
it) w.r.t. the considered group will then yield a symmet- 
ric h which still constitutes a parent Hamiltonian. 

III. LOCALLY SYMMETRIC MPS 

In this section we analyze the implications of a given 
symmetry for a MPS. First, we show that the symmetry 
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transfers to the Kraus operators — generalizing the find- 
ings of [U [3] . In a second step we show that the symmetry 
in the Kraus operators imposes that they are essentially 
uniquely defined in terms of Clebsch-Gordan coefficients. 
Finally, for the special case of SU{2) one can simplify 
even further and analyze the qualitative differences be- 
tween integer and semi-integer spin. 

A. Characterization of symmetries 

It was demonstrated in |3] that the Kraus operators 
which describe any injective state symmetric under a 
group G fulfil the condition J2i'^ij^i = ^g^j^g' where 
u and U are representations of G. We provide in this sec- 
tion a generalization in which injectivity is not required. 
The N appearing in the proof must be sufficiently large 
to obtain property [l] after collecting spins. 

We start by proving the case of discrete symmetries, 
extending the demonstration to continuous groups below. 

Theorem 5 (Discrete symmetries). Let {Ai}f^^ be the 
Kraus operators which describe a locally invariant MPS 
lip) with respect to a single unitary u, i.e. u^^\tl;) — 
e'^lV')- Then, the symmetry in the physical level can 
be replaced by a local transformation in the virtual level. 
This means that there exists a unitary U - which can be 
taken block diagonal with the same block structure as the 
A 's in the MPS and composed with a permutation matrix 
among blocks, i.e. U = P{(BbVb) - such that 

^u,jA^ = WUA,U'' (3) 

3 

With W = ebC^^-'lfc. 

Proof. We follow here a reasoning as in the proof of ^ 
Lemma IV. 4]. We collect the spins in five different blocks, 
each one of them with property [l] Applying u®'^ gives 
us the same MPS (we incorporate the global phase in the 
new matrices) with different matrices S's, but with the 
same block diagonal form and also (after gathering) with 
property [l] We now require the following lemma, which 
is demonstrated below. 

Lemma 6. For each block in the A 's, for instance the 
one given by matrices Aj, there is a block in the B's, 
given by matrices B}, which expands the same MPS. 

Since both are now canonical forms of the same in- 
jective MPS, by O Theorem 3.11], [l^, they must be 
related by a unitary and a phase: ViA\vl — e^^^B}, 
which finishes the proof of the theorem. 

Let us prove now the lemma. By using property [T] and 
summing with appropriate coefficients, it is possible to 
show that there exists a block diagonal D x D matrix 
X ^0 such that 

tr [Al ■ ■ ■ Al] = tr [XB,, ■ ■ ■ , Via, ■ • ■ , is 



Since X ^ 0, there exists one block, let us say Xi, 
different from 0. Then, summing with appropriate coef- 
ficients again we get that there exists a matrix y 7^ 
such that 

tr [YAlAlAl] =tr [X^BIbIb}^ , Vi3,Z4,i5 

We can now argue as in ^ Lemma IV. 4] to conclude the 
proof. □ 

If we have now a symmetry given by a compact con- 
nected Lie group G, that is, ([ij holds for any g £ G and 
a representation g ^-^ Ug, we obtain the following. 

Theorem 7 (Continuous symmetries). The map g 1^ Pg 

is a representation of G and therefore the trivial one. The 
maps g i— > e'^s and g 1— > Vg are also representations of 
G. 

Proof. Let us start with the map g ^ Pg. From eq. ^ 
we get 

'^g2giUg2giAhUl^g^ = uff^^ Aj = 

' (4) 
Y.^'fA^^ = '^g^'^g^.P^.^g^Vg^AnVlvl 
]k 

where Wg-^^p^^ is the same unitary as Wg-^ but with the 
blocks permuted according to the permutation Pg^ . Since 
Pg'Wg — Wg^p^,Pg' s^ud Wg commutcs with all other 
terms appearing in eq. Q, we can multiply successively 
and use property [l] (with L the required block size) , to 
get, for all TT, > L and all X block-diagonal, 

K.grUg.g.^Ulg, = (Wg^Wg^.pJ-Ug^Ug.XUlUl . 

(5) 

By taking X = If, for each block b, we get that Pga^si 
must be Pg^g^. But since we are assuming the group G 
connected, this in turn implies that Pg = 1 for all g. 
With this we can split equation ^ into blocks to get, for 
each h, each n > L and each matrix X, 

"^3251^^ "^3231 ^ * g2^ gi^^ ^ gi ^ g2 v"; 

Taking X = 1 we obtain 

In particular, when n = L, we get that 
L{e''g^gJ = L(6'^^ +OgJ + 2fco7r and when n = L+l that 

+ ^)Kg^) = + + + 2^1'^- Gathering 

both results, the L can be removed and we obtain 

= <+^92+2(fci-fco)^. 

Finally, to show that g 1-^ V^*" is a representation, it is 

enough to notice that eq. |6j) implies that ^g^'i^^g^K^^'sgi 
commutes with every matrix. □ 
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FIG. 2: The unitary Ug applied on the physical level is reflected 
in the virtual level as a pair of unitaries Ug . 



A trivial consequence of these theorems is the fact that 
having an irreducible representation Ug in the virtual 
level implies that the MPS has to be injective. We give 
an alternative proof of this fact in the appendix without 
having to rely on the MPS canonical form. There we 
analyze also when the reverse implication holds. 

B. Uniqueness of the construction method 

Once the theorem which provides the condition that 
the Kraus operators must fulfil in order to generate in- 
variant MPS has been established, the next step is to 
prove that they can always be constructed by means of 
Clebsch-Gordan coefficients. To do that, it is more con- 
venient to work with the map V defined in (|2]). From 
the definition it is clear that the condition ) TiU^^Ai = 
UgAjU^ reads then Ug (E) UgV = VUg. Notice that we 
have removed the dependence on the phase. By Theo- 
rem |7] this can be done for groups with a complex enough 
structure, as SU{2), for which there is no non-trivial one- 
dimensional representation. 

Given a compact group G, the tensor product of two 
irreps -we are choosing a single representative for each 
class of equivalent irreps- can always be decomposed as 
a direct sum of irreps 

Ug(EVgC ^C^Cl 

i 

where C is a unitary whose elements are called Clebsch- 
Gordan coefficients. In what follows we will denote by 
: C''' (g) the matrix associated to the re- 

striction of C to the d^-dimensional invariant subspace 
Hi associated to the irrep c^, with d, d' being the dimen- 
sions of the representations Ug and Vg respectively. 
We are interested in possible solutions of 

Ug (8) Vgfl — flWg . (7) 

where Ug,Vg,Wg are irreps of a given compact group G. 
It is clear that taking 

= (8) 

i 



does the job if we sum over i's corresponding to equiva- 
lent representations = Wg. The next lemma guaran- 
tees that this is all. 

Lemma 8. All possible solutions of Equation ^ are 
given by 

Proof. Any verifying eq. (|7| gives 
n'^n = Wg9)nwl 

which means by Schur's lemma that f2^^f2 = al and we 
may assume that, if there is a non-zero solution, it can 
be taken an isometry. Moreover, introducing V = C^fi, 
which verifies V^V = 1, one has 

vwg = V (9) 

From there one gets that P — VV^ is a rank d pro- 
jector {d the dimension of the representation Wg) that 
commutes with (©iC^) for all g. By Schur's lemma, it 
is supported on ®iHi with i's such that c* = Wg and in 
this subspace it is of the form 

Pi02ld |/32pld ••• - |/?>(/9|®ld . 

This implies that V — \(3) iS)W for a given d x d unitary 
W. But if we substitute this in (|9|, since we are assuming 
a unique fixed representative for each class of equivalent 
representations, we get W = Id and = J^i Pi'Pi- ^ 

From this we can now conclude: 

Theorem 9. Let us consider a group G and two repre- 
sentations Ug (irrep) and Ug — [/^' . Then, the struc- 
ture of all possible maps V fulfilling Ug ® UgV ~ VUg is 

. ' . ' ^ . ' (10) 

KctniVj^^ an2V^^ ■ ■ ■ a„„y^"/ 
where V^^ is a solution, according to Lemma^ to U^^ (8) 

C. The case of SU(2) 

Let us apply the results of the previous section to the 
case in which G — SU(2). Our construction is a natural 
generalization of the one used in [Tl [TB]. 

We consider from now on irreducible representations 
Ug of the symmetry on the physical spin. Nevertheless, 
a substantial part of the results can be straightforwardly 
extended to the reducible case. Hence, we are interested 
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in analyzing the restrictions that SU(2) impose in the 
general solution given by Theorem |9] to the equation 



{U(^J)V = vu 



(11) 



where, with some abuse of notation, J is the SU{2) ir- 
rep corresponding to spin J and U = (ii ® . . . ® z„ ® 
si . . . ® Sm) is the virtual representation composed of 
n integer irreps and m semi-integer irreps. Note that in 
the Clebsch-Gordan decomposition of SU{2) all represen- 
tations appear with multiplicity one. Therefore there is 
only one term in the sum in (|8| . At this point one should 
distinguish the cases of J integer or semi-integer. If J is 
integer, zero is the only solution to {ij ® J)fl = fisj. and 
(sfc (g) J)ri — riij for all j,k, and we get in (10) a block 
diagonal structure: 



V 





V 







a"!/'" 











n+1 ^si 



n-\-m 



VsrJ 



The paradigmatic example in this case is the AKLT 
state [4|, which corresponds to the case oi J — 1,U — \ /2 
in (11). In [1^ , the authors generalized the AKLT model 



to arbitrary integer J and U irreducible. We will call 
the resulting MPS FNW states. It is shown in [T how 
for U = ^ FNW states are unique ground states of frus- 
tration free nearest-neighbor interactions. An alternative 
construction focused on the restrictions imposed by the 
SU(2) symmetry on the density matrix instead of the 
Kraus operators can be found in ^TE\ . 

If J is semi-integer, zero is the only solution to {sj (E) 
J)fi = and {ik ® J)^ = ^ij for all j, fc, and we get 
in ([To|) an off-diagonal structure: 



V = 
/ 



■^n+l " si 



n+1 











n-|-m 



ysm\ 







It is clear that the virtual representations must be re- 
ducible now, which is very much related to the Lieb- 
Schultz-Mattis theorem, as we will show in the following 
section. The paradigmatic example in this case is the 
Majumdar- Ghosh model [5j, which corresponds to J = ^ 
and ?7 = 1 0. A generalization of this model for the 
case of arbitrary J and U — F0O, was recently proposed 
in M- 

In general, it is possible to find a set of representa- 
tions which fits into any model with SU(2) symmetry, 
for instance [T71 [HI |^ . 



IV. LIEB-SCHULTZ-MATTIS THEOREM 

The Lieb-Schultz-Mattis theorem states that, for semi- 
integer spin, a S'?7(2)-invariant ID Hamiltonian cannot 
have a uniform (independent of the size of the system) 
energy gap above a unique ground state. That is, symme- 
try imposes strong restrictions on the possible behaviors 
of a system. In this section we want to go a step further 
and analyze which implications one can obtain from hav- 
ing a single symmetric state in a semi-integer spin chain. 
By restricting our attention to the class of MPS we will 
show 

Theorem 10. Any MPS with an SU{2) symmetry in 
the sense of |Ip with Ug irrep and even physical dimen- 
sion d cannot be injective. By Theorem 11 of ^2! this 
implies that it cannot be the unique ground state of a,ny 
frustration free Hamiltonian. 

Proof. Let us assume that the MPS is injective and prove 
the theorem by contradiction. Theorems [S] and [7] guar- 
antee that 



A, = UgAkUl . 



(12) 



We consider 

k — 1, . . . ,d. Then, eq. (|12| gives 



- e'J' with (J,),-fc = 5,, fe(fc - (rf+ l)/2). 



2* "P^Ak = UAkU'' 



(13) 



for a unitary U and (pk half-integer. We finish by proving 
that if N is odd, tr{Aki ■ ■ ■ A^^) = and hence the MPS 
cannot be injective. From ( 13 1 we get tr( • • • Ak„ ) — 
unless J2i=ifki — N{d+ l)/2. The latter is, however, 
impossible for N odd as then the l.h.s. is integer whereas 
the r.h.s. is half-integer. □ 

From the proof one may get the impression that 
only U{1) symmetry is required, and this is indeed the 
case if the generator of such symmetry has eigenvalues 
— r7i/2, m/2 as above. The next example shows that 
this is, however, not true for any U{1) symmetry, which 
in turn shows that a larger symmetry like SU{2) is re- 
quired for the Lieb-Schultz-Mattis theorem. 

Example 11. Let us consider a local symmetry gen- 
erated by G = e^^^ for a hermitian matrix H. Let 
us choose the physical dimension d = — D, which 
is always even, and the set of Kraus operators JC = 
{A(i.j) = j|, i^j}. Select ai,...,aD G M such 
that — Qfj 7^ if i ^ j and H the diagonal ma- 
trix H — ~ (which has in addi- 
tion only non-zero eigenvalues). With Ujs = e^^^ where 
Q. = diag[ai . . .ao] it is clear that 



gi/3(ai-Qj)^ 



SO the MPS generated by means of the Kraus operators K. 
has the local symmetry G. Moreover, the MPS is trivially 
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injective when D > 3. We can prove this by choosing 
arbitrary k and k'. Since _D > 3, we can always find an 
I such that k' ^ k and then \k){k'\ = \k){l\l){k'\ = 

Let us remark that this counter-example is applicable 
to spin > |. Indeed, one can prove Theorem for U{\) 
and spin \, which is the content of the following propo- 
sition. The case of spin 3/2 remains an open question. 

Proposition 12. // |(f>) is an MPS with physical dimen- 
sion d~2 and invariant under U{1), then |<i>) cannot he 
injective. 

Proof. We will show it by contradiction. By choosing a 
basis where the physical unitary u is diagonal, the con- 
dition on the Kraus operators becomes 

where H is the hermitian generator of the symmetry. Let 
us expand the expression for infinitesimal angles 

which is the equation of eigenvalues for the operator 
L(») = [i?, •]. This can be transformed into an ordinary 
eigenvalue equation for the matrix operator L = H ® 
1 — 1 ® H . The diagonalization can be easily performed 
by taking the spectral decomposition of H = '^ifJ-iPi, 
where Pi are orthogonal projectors. It straightforwardly 
follows that the eigenvalues of L are Ay = I-M — /ij and 
the corresponding eigenoperators fulfil Aij = PiAijPj. 

Let us focus now on the case d = 2. Then, we have 
that Ai = PiAiPa and A2 = PpA2Pj for some a,f3,j. 
If /3 = 1 PiX = X for all X e span{Aj, ■■■A,J and the 
MPS cannot be injective. The same happens if a = 7. 
So let us assume that /3 ^ 1 and 7 7^ a. Now if a 1, 
we have Ai = PiAiPi, A2 = (1 - Pi)A2{l - Pi) and the 
MPS is block diagonal and hence non-injective. The same 
happens if /3 = 7. So a ^ 1 and /? 7 and this gives 
A\ = Q = A\ which implies that spanj^i^ ■ ■ ■ Ai^} = 
span{AiA2AiA2 • • • , A2A1A2A1 • • • } has dimension < 2. 

□ 



V. GENERAL CONSTRUCTION OF 5(7(2) 
TWO-BODY HAMILTONIANS WITH MPS 
EIGENSTATES 

We have seen in Definition |3] a way, called the par- 
ent Hamiltonian method, to construct local SU{2)- 
symmetric Hamiltonians with MPS as eigenstates. In 
this section we first prove that this method is the most 
general one to find Hamiltonians having a given MPS as 
local eigenstate, that is, being an eigenstate of each local 
term in the Hamiltonian. Then, we show examples (in- 
cluding the AKLT and Majumdar-Ghosh states) of MPS 
that are excited eigenstates of local two-body transla- 
tionally invariant S'?7(2)-symmetric Hamiltonians. More 
examples are then provided in the appendix. 



A. Completeness of the parent Hamiltonian 
method 

Theorem 13. Given an MPS \ip), any translational in- 
variant Hamiltonian having it as a local eigenstate is of 
the form al -\- H where H is a parent Hamiltonian for 
in the sense of Definition^ 

Proof. Let us call h the local hamiltonian. By hypothesis 
of local eigenstate, 

hp = Xp (14) 

for certain A G R. This implies [p, h] = and hence one 
can find a set of projectors V = {Pi, i = 1, . . . ,r\ ^ ■ = 
1} such that we can decompose both p and h by means 
of them, i.e. h = J2i O'iPi and p — X^jgc ^i^i^ where C 
represents the set of projectors which describe the sup- 
port of p. Using eq. ( [14^ with this decomposition gives 
that fli = A for alH G C and hence 

iec^ tec 
= X] - ^)P^ + ^1 ■ 

Then, the translational invariance hamiltonian is H — 
'}2,jT^{h) (g) licst, where r is the translation operator. 
The theorem follows from replacing the result for the 
local hamiltonian and comparing this with Definition [3] 
of parent hamiltonian. □ 

This Theorem shows that, given an MPS j?/;), looking 
for all possible parent Hamiltonians of interaction length 
k is equivalent to look for all possible solutions to the 
equation 

V^^^Ap^'^), (15) 

with A = tr [/ip''^''] . The next lemma gives yet another 
equivalent formulation, which is the one we will use in 
the sequel. 

Lemma 14. Given a Hermitian matrix h and a density 
matrix p, hp — Xp if and only if 

tr [h^ p] -tr [hpf ^0 . (16) 

Proof. One implication is clear. For the other, let us 
write (h) for tr [hp] 1 . By assumption 

tr [{h - {h))^p] = tr [h'^p] - tr [hpf ^ 0. 

So p^l'^{h — {h))^p^/^ — 0, since it is a positive operator 
with trace 0. This implies that {h — {h))p — and hence 
hp = Xp. □ 

With this at hand we can systematically search for 
MPS that are excited local eigenstates of SU{2) invari- 
ant Hamiltonians with two-body interactions. We will 



proceed as follows. We start with a given SU{2) sym- 
metric MPS lip) and fix the interaction length n. Then 
we look for possible solutions to Eq. ( 16 1 of the form 



2J 

i<j<n a— 1 



(17) 



to ensure SU(2) symmetry and two body interactions in 
the Hamiltonian. Finally, to guarantee that the MPS 
is an excited state, we will find another SU(2) symmetric 
MPS with less energy that will act as a witness. In the 
next section we will illustrate this procedure starting with 
!■(/;) the AKLT, the Majumdar-Ghosh state, and general- 
izations. Throughout we work in the thermodynamical 
limit N oo. 



B. Examples of SU(2) two-body Hamiltonians 




1. Spin 1 

Let us consider the AKLT state as a first example. Its 
Kraus operators are A_i = — \/2cr~, Aq — cr^, Ai — 
V2a+. 

In the case n = 2 the only solution to Eq. (16 1 is the 
AKLT Hamiltonian. In the case n = 3, the solutions are 
given by 

h = (-3wi + V2 + 3v3){Si o ^2) + v^^Si o 82)'^+ 
^(-3wi -I- V2){Si o S3) - ^(-3ui + V2){Si o 83)'^+ 

V2{S2 O S3) + Vi{S2 O 83)^ 

where the eigenvalue corresponding to the AKLT state is 
7vi — Sv2 — 2^3. The total translational invariant Hamil- 
tonian is then 



H = ^(-3«i + 2v2 + Sv3)iS, o S,+i)- 



{Vl + V3){Si O S'i+i)^ -I- 2^-3'"^ + V2){Si O Si+2)- 

^(-3ui -I- v2){Si o 81+2)'^ 

which contains the usual AKLT model. It is not difficult 
to check that there is a region in the parameter space 
where the AKLT state is still the ground state of this 
Hamiltonian. To find regions where it is an excited eigen- 
state we will use as a witness the 8U{2) symmetric MPS 
associated to the virtual representation |® 5 (see Section 



nil. The result is plotted in fig. |3j where one sees the 
existence of points in this family of spin 1 Hamiltonians 
for which the AKLT state is an excited state. 

Note that it is possible to perform a change of variables 
in the total Hamiltonian, for instance a |(— 3wi + 
V2) and b ^ Vi + V3, such that it depends only on two 
parameters. However, the number of parameters that 



FIG. 3: Space of parameters of the local Hamiltonian h for the 
AKLT state and n = 3. The orange volume represents the points 
where the state is the local (and hence the global) ground state. 
The green volume represents points corresponding to excited states 
detected with the witness | © i . 



^\ j 


1 


1 


3 


2 


5 


3 


7 


J 


2 


2 




2 




2 


1 


2 


1 












2 




5 


3 


■ 


■ 


■ 


■ 


3 






4 


2 


2 


1 


■ 



TABLE I: Table of results for FNW states with physical spin J 
and virtual spin j. The numbers in the table are the number of 
parameters the obtained families of Hamiltonians depend on. The 
I represent the cases for which no solution was found. 



the local Hamiltonian h depends on cannot be reduced, 
which means that there are non-physical parameters in it. 
In Fig. |4] we have represented the problem above (n = 3 
and AKLT state) in terms of the physical parameters. 
The positive axis b corresponds there to the usual AKLT 
Hamiltonian. 

Concerning FNW states, that is integer spin J and 
virtual irrep j, we have performed an exhaustive search 
and table |l] gathers the main results. The study has been 
carried out by increasing n and studying the number of 
parameters which the family of Hamiltonians depends 
on (notice that the case of interaction length n contains 
the case of interaction length n — 1). We have increased 
n until the number of parameters stops growing. In all 
the cases considered in the table, a saturation occurs 
when n > 3, i.e. considering more than 3 particles does 
apparently not add new Hamiltonians. 

Let us also introduce a new state of spin 1 with virtual 
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FIG. 4: Space of physical parameters of the global Hamiltonian H 
corresponding to n = 3 and the AKLT state. The orange points 
represent where the state is the local (and hence the global) GS. 
The green surface represents points corresponding to excited states 
detected by means of the witness | © j . 



spin 1, given by the Kraus operators 
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The total translational invariant hamiltonian which has 
this state as eigenstate is 



This state is injective and a local excited state. The fact 
that this state is an excited state of the global hamilto- 
nian can be checked as above by means of the witness 
1® 0. 



2. Spin \ 

Let us consider now the the Majumdar-Ghosh state as 
an example with semi-integer spin. The Kraus operators 
are now 



A_i 




Ai 





1 

75 0. 



As in the previous case, we do not find any solution for 
n = 2 and only the Majumdar-Ghosh Hamiltonian for 
the cases n — Z and n — A. For n = 5 the solutions to 



Eq. ( 16 1 are given by 



h = {vi — V2 + V4){Si O S2) + (Vi — V2+ V4){Si o 83) + 
V3{Si O §4) + V3{Si O 5^5) -I- V4{S2 O 5^3) 
+ {-Vl +V2 + V3){S2 O 5*4) + V3{S2 O + V2{S3 O 5*4)-^ 

Vi{S3 0S5) + Vi{S4 0S5) (18) 

and the energy associated to the state is — |(tii -I- U4). 
The total Hamiltonian H — J^i Ti{h) is given by 

H = 2(t;i -I- V4){Si o 5.1+1) + (^^1 + ^^3 + V4){Si o 81+2) 

i 

+ 2v3{S,oS,^3) + V3{S,oS,+4) (19) 

As in the AKLT case, by means of a change of variables 
a ^ V3 and b —* vi + V4, the number of physical param- 
eters in the total Hamiltonian is 2, compared with the 
four parameters the local Hamiltonian depends on. The 
Majumdar-Ghosh state is an excited local eigenstate for 
a region in the space of parameters, which in this case is 
detected by the witness | © 1 ffi 0, as shown in fig. [s] The 
usual Majumdar-Ghosh Hamiltonian ^24j corresponds to 
the positive axis b. 

3. Spin I 



Let us consider as final example the SU (2) symmetric 
MPS corresponding to spin | and virtual representation 
I 0. For n = 3, the solutions to Eq. ( 16 ) are given by 



h — V3{Si o 52) + V2{Si o 52)^ + vi{Si o 52)^4- 
{2vi -V2 + V3){Si o 5*3) -I- (4i;i - V2){Si o ^3)^ 

-I- Vi{Si O 53)^ + V3{S2 O S3) + V2{S2 o S3Y + 

vi{S2oS3f (20) 

and the energy associated to the MPS is in this case 
— g|(165'i;i — 60u2 -t- 161)3). The global Hamiltonian reads 
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FIG. 5: Space of physical parameters of the total Hamiltonian for 
n = 5 associated to the Majumdar-Ghosh state. The orange points 
represent where the state is the local (and hence the global) ground 
state. The green surface represents points corresponding to excited 
states detected by means of the witness g © 1 © 0. 



now 

i 

+ 2vi{S, o S,+if + (2wi --V2 + V3){S, o S,+2)+ 

(4i-i - V2){S, o S,+2f + vi{S, o S,+2f (21) 

It is remarkable that in this case there are no spurious 
parameters in the local Hamiltonian h. Considering the 
family of states whose virtual representation is | © 1 © 
as a witness, it is possible to demonstrate that there is 
a region in the space of parameters of the Hamiltonian 
for which the MPS is an excited eigenstate, as shown in 
Fig.|6] 

VI. CONCLUSIONS 

Despite the fact that all our results are restricted to 
the family of TIMPS, their relevance is manifested by the 
fact that those states approximate all ground states of 1- 
dimensional Hamiltonians with short range interactions. 
Thus, one would expect that the properties derived for 
MPS would be relevant in a more general context. More- 
over, due to their simplicity, MPS can be then thought 



as a 'laboratory' where to search for some generic math- 
ematical and physical properties of states that are rel- 
evant in 1-dimensional spin chains. Later on, one may 
use more powerful mathematical methods to try to ex- 
trapolate those properties to general spin chains. Fur- 
thermore, many of the techniques used in the present 




FIG. 6: Space of parameters of the spin | model. The orange 
points are obtained numerically and they represent values of the 
parameters where the MPS state is the GS. The green volume rep- 
resents points corresponding to excited states detected with the wit- 
ness I © 1 © 0. 



work are amenable of an extension to higher spatial di- 
mensions, where PEPS play the role of MPS. In Ref. [5] 
some first results in this direction were derived, which 
will be generalized in a further publication. 
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APPENDIX A: RELATIONS BETWEEN 
IRREDUCIBILITY AND INJECTIVITY 

In this appendix we give a direct proof of the fact that 
an irreducible representation in the virtual level of a sym- 
metric MPS implies that the MPS is injective. We also 
see that the reverse inclusion is not true in general, but 
it holds under some conditions on the Kraus operators. 

We have to recall that, given a set of Kraus operators 
defining an MPS K, = {Ai, . . . , Aj,}, we can define an as- 
sociated completely positive map K{X) = '^i^iAiXAj. 
The symmetry in the MPS transfers then to the covari- 
ance of the channel, that is, E{UgXWg) = UgE{X)Ul for 
all X. It is shown in [Tl that if E is trace preserving 
and has 1 as its unique fixed point, then the MPS is in- 
jective. Moreover, it is trivial to see that if E is the ideal 
channel {E{X) = X for all X), then the MPS is a prod- 



uct state. Therefore, the desired result that irrep implies 
injectivity is a consequence of the following theorem. 

Theorem 15. Let us take a completely positive map 
E : A4d — > -Md that is covariant for an irrep of a 
compact connected Lie group G. Then, either E is the 
ideal channel or it is trace preserving and the identity its 
unique fixed point. 

Proof. Let us consider a fixed point A of E. Then UgAU^ 
is also a fixed point because of the covariance. Therefore, 
integrating under the Haar measure and using Schur's 
lemma, 1 is also a fixed point. A similar argument shows 
that E is also trace preserving. 

Now we can apply Liiders' theorem [7], which ensures 
that the set of fixed points of E coincides with the 
commutant JC' of the set of Kraus operators of E. This is 
trivially a C*-subalgebra of Md- Moreover, we know by 
the classification of the C*-subalgebras in Md that there 
exists a unitary V £ Md such that VVV^ = ©i(M„^ ® 

The equivalent representation Vg = VUgV^ is also an 
irrep and fulfils that VgAV^ — A. This means that the 
block structure of A remains invariant under the action 
of Vg by conjugation. Now we use that 

VgAV^ C A^[J,A]C A for all generators J. (Al) 



This implies that J has the same block structure as A. 
If there is more than one block, the representation is 
reducible. If ^ = M„ ® 1„/, then we use again eq. ( Al I: 



The Schmidt decomposition allows us to take J 
Ai^Bi where the Bi 's form a basis of M„/ , with Bi 
1. Then, 



eq. 



AT] gives that I]J^j, M„] (g) Bi = C (g) 1, 
which implies that Ai is proportional to 1 for all i > 2. 

1 and hence Vg = Vf (g V2 , 
1 or ^ = Mjv (which im- 

□ 



This gives J =1<^X + Y ( 
which is reducible unless A 
plies that E is the ideal channel). 



Although the implication in the opposite direction 
could also seem true, it is not, as shown by the following 
example. 

Example 16. Let us consider the family of SU (2) sym- 
metric MPS of spin 1 with a reducible virtual represen- 
tation 5 ® I given by the following maps (see Section 



III) 



V 



( e*"ii cos 6*1 V? e*"" sin 6*2 \ 

2 2 

. e*"" sin OiVi e*"^^ cos 6*2^1 ; 

\ 2 2 ' 



It is not difficult to check that the MPS is injective except 
in particular directions in space, such as those for which 
the isometry breaks into blocks, i.e. 6i = n^. 

Although the equivalence is not true in general, we can 
still give a sufficient condition which applies, for instance, 
to the AKLT and other FNW states. Let us recall from 
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[3] or Theorems [5] and [7] that an injective symmetric MPS 
verifies 

Y,uf^A, = e^'^U,A,Ul, (A2) 

i 

where in addition one may ask for ^l^i = 1 

Proposition 17. // Ug is irreducible and {AlAj}ij 
spans the whole space of matrices, then the virtual repre- 



j = l,e=l: 



sentation Ug of ( A2) is also irreducible 



Proof. From (|A2|) one gets 



E 



-9 9 



^n^«2 — UgA^^^Aj^U^ 



Integrating now with respect to the Haar measure, the Ihs 
is simphfied by the irreducibihty of Ug and the orthog- 
onality relations. The result is Sj^j^J^iM^i — ^jij2'^- 
This means that J^UgXU^ (x 1, VX S A4d, since we 
can span the complete space of matrices. But this im- 
plies that Ug is an irrep by means of the inverse of Schur's 
lemma. □ 



H = '^{Si o - (Si o Si+2) - {Si o 5^+2)^ 



No solutions found (with n < 4) for j |, 2, |, 3. 



3. Spin J = f 

• j= |©0, e=-i|(165wi-60v2 + 16w3): 

i 

2vi{Si o S'l+i)^ + {2vi ~V2 + vz){Si o 5^+2)+ 

(4di - V2){Si O Si+2Y + O Si+2)^ 



APPENDIX B: LIST OF PARENT 
HAMILTONIANS 

The following lists SU (2)-invariant two-body Hamilto- 
nians for which the MPS with physical spin J (irrep) and 
virtual spin j is an exact eigenstate with energy e. 



1, e: 



495 . 

64 ■ 



1. Spin J = i 



No solutions found (with n < 4) for j = | ® 1,| 



1, ^ 2, ^ 2. 



H = Y2{vi+V4){SioSi+i) + {vi + V3 + V4){SioS^+2) 

i 

+ 2vi{Si o S^j+a) -I- V'i{Si o Si+i) 
• No solutions found (with n. < 6) for j ~ 5 1; 



2. Spin ,7 = 1 

• j = 5, e = Twi - 3v2 - 2V3: 

H = Y^{-ivi + 2v2 + 3v3)(5, o 

i 

{vi + V'i){Si o Si+i)'^ + ^(-3^;l -I- V2){S.i o Si+2)- 

1 ^ ^ „ 

-(-3ui -I- V2){Si O S'j+2)^ 



4. Spin J = 2 



j = 1, e = (-6986wi-|-778w2-62w3-|-1260u4-90w5): 



iJ = Y{2mvi-63v2+2Av3-792vi+63v5){S,oS,+i)- 

i 

(133W1 - 14w2 + 2w3 - 133i;4 + 14w5)(S, o Si+if+ 
{v2 + v^){Si o Si+i)^ + {vi + Vi){Si o Si+i)'^+ 

i^-vi - 91u2 + y W3)(jj o <Si+2)+ 
5719 301 43 ^ ^ 2^ 

, 665 35 5 ..-^ .3 

( YS^I + y«2 - YQ^3)(Sr O S^+2) + 

(-^^'1 + 1^2 - Y^Vs){Si O 8^+2^ 



H = ^(580vi)-80t;2+10t;3-330?;4+30w5)(<5iO<5i+i)+ 

i 

(91wi - nV22vj. - QlVillV5){S^ O S.,+ if + 

{v2 + v5){St o Si+if + (t'l + V4){S^ o 5^+1)^^+ 
^(2275t;i - 275V2 + 25v3){Si o 5^+2)+ 

^(455z;i - 55v2 + 5vs)iSi o 5^+2)'+ 

4(-455wi + 55v2 - 5v3){S^ o 8^+2)^+ 
10 

• No solutions found (with n < 4) for j = 2, |. 
5. Spin J = 3 



Solutions (mostly cumbersome ones) were found for 
= l(n = 3), j = 2{n = 2) and j = 5/2(n = 2). 



